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Abstract. The interstitial oxygen centers in silicon and germanium are reconsidered and 
compared in an analysis based on the first-principles total-energy determination of the 
potential-energy surface of the centers, and a calculation of their respective low energy exci- 
tations and infrared absorption spectra. The total-energy calculations reveal unambiguously 
that interstitial oxygen is quantum delocalized, the derealization being essentially different 
in silicon and in germanium. Oxygen in silicon lies at the bond center site in a highly an- 
harmonic potential well, whereas in germanium it is found to rotate almost freely around the 
original Ge-Ge bond it breaks. This different derealization is the origin of the important 
differences in the low energy excitation spectra: there is a clear decoupling in rotation and 
vibration excitations in germanium, giving different energy scales (1 cm -1 for the rotation, 
200 cm -1 for the V2 mode), whereas both motions are non-trivially mixed in silicon, in a com- 
mon energy scale of around 30 cm -1 . The calculation of the vibrational spectra of the defect 
reveals the existence of vibrational modes (related to the V\ mode) never been experimentally 
observed due to their weak infrared activity. It is found that the combination of these modes 
with the well established 1/3 asymmetric stretching ones is the origin of the experimentally 
well characterized modes at frequencies above the 1/3 mode frequency. 

I. Introduction 

It is broadly accepted that interstitial oxygen in silicon (Si:Oi) [1] and in germanium (Ge:Oi) 
[2-5] represent very similar centers in their geometry and dynamics. The aim of this paper is 
to show the contrary: the quantitative differences in the low-energy experimental results for 
Si:Oi [6] and Ge:Oi [3] are shown to correspond to fundamentally different realities, that, in 
turn, lead us to predict substantial differences in their infrared absorption spectra. The basic 
physics of both centers as well as their similarities and differences are analyzed and discussed 
in this work, the essence being graphically depicted in Fig. 1. Oxygen breaks a covalent bond 



Fig. 1: Schematic representation of the 
oxygen derealization in Si:Oi and Ge:Oi, 
in the plane perpendicular to the Si-Si 
(Ge-Ge) axis. 



between two semiconductor atoms and remains quantum delocalized between them. The 
derealization is essentially different in Si:Oi and in Ge:Oi, both being again different from 
the trivial harmonic motion around a well defined geometry. In the case of silicon [7,8], 
oxygen is delocalized in the close neighborhood of the original silicon-silicon bond center 
(BC) position, whereas in germanium the defect can be very precisely accounted for by an 



almost free rotor. As a consequence of the different geometries the low energy excitation 
spectra are quite different. In addition, in the infrared spectral region, interstitial oxygen 
in silicon induces a backbond-stretching resonance at 517 cm -1 which has no counterpart in 
Ge:Oi, while in germanium there is a z/ 2 bending mode in the infrared that in Si:0^ appears 
among the low-energy excitations. 

In this work we address at these differences between Si:Oi and Ge:Oi from different and 
complementary points of view. With the help of total energy calculations, we determine the 
"equilibrium" configuration of Si:Oi and Ge:Oi and analyze the low energy excitations experi- 
ments. The atomic vibration and infrared spectra associated to the defect are also calculated. 
We find vibrational modes which have never been directly observed experimentally. We show 
that these modes, in combination with the high frequency stretching modes, are responsible 
for the experimentally observed modes induced by oxygen at 1751 cm -1 and 1260 cm -1 in 
Si:Oi [9] and Ge:Oi [2], respectively. 

II. Atomic configuration and far infrared analysis 

In Fig. 2 the results of the analysis of the low energy excitations of the interstitial centers are 
shown. In the case of Si:0; the theoretical results of Yamada-Kaneta et al. [7] along with the 
experimental results of the far infrared measurements of Bosomworth et al. [6] are displayed, 
supporting the image of Fig. 1. In the case of Ge:Oi, the experimental results of Gienger 
et al. [3] are drawn along with the results of our hindered rotor theoretical analysis, where 
the Hamiltonian of an elastic rotor [10] (H = Bl 2 — Dl 4 , I being the angular momentum) is 
perturbed by an angular-dependent potential H' = (A/2) cos(6(/>) with an amplitude A m 0.6 
meV. The model is solved exactly to the required precision, giving a remarkable agreement 
with experiments, which supports the picture of Fig 1. It has to be stressed, however, that 
higher harmonics of the angular potential affect very little the displayed results, leaving the 
hindering potential (and thus the angular energy barrier, and the position and number of 
minima) indetermined in that respect (the measurement of the splitting of the / = ±6 levels 
would resolve the indetermination) . This is illustrated in Fig. 2, where an added cos(12</>) 
term is considered. Note the possible qualitative differences. 



Fig. 2: Low-energy excitations of 
Si:Oj and GeiO^. For Si:0^ far-infrared- 
absorption data [6] are compared with the 
fit of Yamada-Kaneta et al. [7] based on 
the anharmonic potential well coupled to 
the Vs vibration mode. Also displayed is 
the renormalized potential for the ground 
state of z/3. For Ge:Oi, the results of 
phonon spectroscopy of Gienger et al. [3] 
are compared with the results of the hin- 
dered rotor model. Also shown are three 
potentials compatible with the data, dif- 
fering in the amplitude of the cos(12(/>) 
term (zero, negative, or positive, respec- 
tively). 



It is important to note in Fig. 2 the different energy scales involved. The angular de- 
localization of oxygen in germanium gives rise to low lying rotation excitations which are 
essentially decoupled from the vibrational modes at a different energy scale (see below). In 
silicon, radial and angular excitations are non-trivially mixed in a common energy scale. 

First-principles total-energy calculations have been performed in the manner described 
elsewhere [8], using a cluster-Hartree-Fock approximation. The results of the calculations, 
carefully considering relaxations up to second nearest neighbor atoms of oxygen, are shown 
in Fig. 3 for the radial dependence. The size of the barrier for the radial motion in Ge:Oj 
localizes the oxygen distance to BC around ro = 0.58 A. The amplitude of the cos(6</>) 
harmonic for the angular dependence of the potential has been calculated to be of a few 
tenths of a meV for Ge:Oi, in qualitative agreement with the model described above. These 
results give additional support to the picture of Fig. 1. In the equilibrium geometry the 
Ge-O-Ge (Si-O-Si) angle is 140° (180°) and the Ge-0 (Si-O) distance 1.70A (1.56A). 

Fig. 3: Radial potential for the oxy- 
gen motion for Si:Oj and Ge:Oi ob- 
tained from first principles relaxing up 
to second nearest neighbors of oxygen 
for every point. There is a flat region 
in the Si:Oi case [8] compatible with 
the potential in Fig. 2. 



III. Vibrational modes and infrared absorption 

The infrared absorption spectrum of Ge:Oi is calculated for a fixed polar angle, in absolute 
units [11] and the results are shown in Fig. 4 (the results for Si:Oi have been published 
elsewhere [8]). Four distinct features appear at 877, 416, 230, and cm -1 , respectively. The 
latter corresponds to the free rotation of oxygen and the other three correspond to the z/3 
stretching, ui bending, and z/ 2 rocking modes of the Ge20 unit [1] (see below). The never ob- 



Fig. 4: Infrared absorption coefficient 
calculated for Ge:Oi for an oxygen con- 
centration [0]= 10 18 cm -3 . The in- 
set shows a magnification of the spec- 
tral region corresponding to the upper Ge 
phonon bands, compared with the equiv- 
alent in Si:Oi. 



served peak appears as a resonance in the germanium continuum and corresponds to the 
radial vibration of the rotor. Indeed, the centrifugal distortion coefficient D = 4B 3 /h 2 u 2 [10] 
in the hindered rotor model gives a value for u> of the magnitude of the frequency of this 
mode (the exact value depends on the amplitude of the cos(120) term). The position of this 
h> 2 feature is not as well defined by our model as the others, being more sensitive to small 
changes in the dynamical matrix. 



Fig. 5: Displacement-displacement 
correlation functions, (uiUj){u) of 
oxygen and neighboring germanium 
atoms in GeiO^. 



Although the asymmetric streching mode at 877 cm -1 has been very well characterized 
experimentally as a main absorption peak at 860 cm -1 , the mode at 416 cm -1 (parallel 
to the 596 cm -1 mode in Si:Oi) has never been observed experimentally. This is due to 
the weakness of this feature and also to the fact that the infrared spectrum of crystalline 
germanium shows a two-phonon infrared absorption in this frequency range [12] that can 
hide the absorption of at the usual concentration. Nevertheless, the combination of this 
mode with the asymmetric stretching at 877 cm -1 has been observed experimentally at 1260 
cm -1 [2], in the same manner as in Si:Oi the non-infrared active mode at 596 cm -1 combines 
with the asymmetric stretching one to give the 1751 cm -1 band [9]. This assignment of 
the combination mode is soundly supported by the agreement of the isotope shifts predicted 
by the theory and the corresponding experimental results. These data have already been 
presented for S'v.Oi [8,13] and will be published elsewhere for Ge:0; [14]. 

The displacement-displacement correlation functions (uiUj)(uj) of atoms in the vicinity 
of the defect are shown in Fig. 5, confirming the assignments made above. The inset of 
Fig. 4 shows the infrared absorption near the bulk continuum for both SkO^ and Ge:Oi. The 
figure shows remarkable differences between these two spectra. The asymmetric resonance 
in Si:Oi at 517 cm -1 , absent in Ge:Oi, is due to the backbonding relaxation caused by the 
incorporation of oxygen. This relaxation is larger in silicon than in germanium. For a better 
understanding of this feature we have calculated the displacement-displacement correlation 
function between the silicon atoms forming the backbonds. Results of the calculations are 
shown in Fig. 6. These results show very explicitly the origin of the modes at 517 and 596 cm -1 
(416 cm -1 in Fig. 5) in silicon (germanium). The mode at 517 cm -1 is a silicon optic mode 
(transverse to the defect axis) , with a higher frequency due to the backbonding compression. 



The mode at 596 cm 1 (416 cm x ) is a (longitudinal) optic mode split-off the band due to 
the presence of the stronger Si-O-Si (Ge-O-Ge) bond in the semiconductor lattice. 

Fig. 6: Displacement-displacement correlation func- 
tions of silicon atoms forming backbonds in Si:Oi. 

IV. Concluding remarks 

Total energy calculation, analysis of the low energy excita- 
tions and calculation of the infrared spectra have allowed us 
to fully understand the similarities and differences of inter- 
stitial oxygen in silicon and germanium. It is found that, in 
spite of the similarity of the defects, there are fundamental 
differences that manifest themselves in the low energy exci- 
tation spectra and appear as subtle but relevant details in 
the infrared absorption. 
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